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ABSTRACT: 

An  existing  computer  simulation  of  a  hydraulic  submarine  steering 
control  system  was  analyzed  to  determine  the  source  of  its  unrealistic 
behavior,  which  included  verifying  mathematical  equations  and  numerical 
parameters,  duplicating  the  simulation  on  another  computer,  and  investi- 
gating the  validity  of  the  assumptions  made  in  its  development.   Rudder 
turning  rate  and  load  pressure  drop  oscillations  which  were  thought  to 
be  unrealistic  were  caused  by  the  neglect  of  significant  rudder  moments 
related  to  the  control  surface  turning  rate.   Knowledge  of  actual  system 
performance  was  found  to  be  incomplete  and  further  research  was  suggested. 
Potentially  beneficial  modifications  in  the  servovalve  and  actuator 
simulators  were  identified.   Continuing  work  in  the  field  of  unsteady 
rudder  behavior  was  suggested  in  order  to  improve  simulator  and  actual 
system  designs. 
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LIST  OF  SYMBOLS 

A  actuator  piston  area 

a  potentiometer  setting 

A  crosshead  area 

CH 

A  area  of  an  incremental  flat  plate  used  in  the  calculation 

^  of  stationary  rudder  damping 

A'  area  of  holes  in  crosshead 

C  loss  coefficient  of  a  short  tube  orifice.   Also  the  chord 
length  of  a  rudder. 

c  radial  clearance  of  a  cross  head  in  its  sleeve 

C  coefficient  of  friction  of  the  stock  bearing.   Also  orifice 
drag  coefficient. 

C,^-,  dimensionless  twisting  moment  coefficient 
Mrl 

C^,-  quasi-stationary  portion  of  C,_, 
MHO  Wn 

D  diameter  of  axial  holes  in  crosshead 

db  dead  band,  fraction  of  maximum 

E  symbol  for  the  Laplace-domain  variable  for  the  error  signal 

e  eccentricity  of  crosshead  in  sleeve.   Also  the  error  signal 

f  function  of  some  yet-to-be-determined  form 

I  polar  moment  of  inertia  of  the  rudder 

K  amplifier  gain  coefficient 

K  flapper  centering  spring  rate 

K_  wand  spring  rate 

K,  pilot  stage  valve  flow  coefficient 

K  main  stage  valve  flow  coefficient 

K,  system  hydraulic  resilience  coefficient 
o 

K_  hydrodynamic  load  torque  coefficient 

K  hydrodynamic  load  torque  coefficient 

o 

K„  lift  coefficient 


K  ^  seal  friction  coefficient 

K  drag  coefficient 

K  mechanical  resilience  coefficient 

L  tiller  length 

1  length  of  various  holes  or  channels 

M  moment 

M^  twisting  moment  generated  by  a  rotating  rudder  in  a  moving 
fluid 

N  number  of  holes  in  the  crosshead 
o 

N  number  of  slots  in  the  crosshead 
s 

P  pressure  in  the  forward,  enclosed  portion  of  the  crosshead 
sleeve 

P^  pressure  in  the  aft,  open  portion  of  the  crosshead  sleeve 

I. 

P  load  pressure  drop 

P,    ,  P^    maximum  load  pressure  drop 

Lraax  Lm 

P  supply  pressure 

Q  hydraulic  fluid  flow  rate 

Q  flow  rate  through  an  annulus  formed  between  two  eccentric 
cylinders 

Q  flow  through  the  rectangular  slots 

Q  flow  through  the  crosshead  holes 

R  rudder  stock  radius 

r  radius  of  the  crosshead  sleeve.   Also  distance  from  the 
pivot  line  of  the  rudder. 

s  Laplace  transform  variable 

S  rudder  span 

T  variable  used  in  the  expression  for  T 


c 


F 


T  equivalent  torque  generated  by  the  crosshead 

CH 

T_         friction  torque 
r 

T  hydrodynamic  torque 

H 


T  pressure  torque 

T  variable  used  in  the  expression  for  T 

s  '  F 

T»  rudder  torque  generated  by  a  rotating  rudder  immersed  in  a 

stagnant  fluid 

T  system  sliding  friction  torque 

U  velocity  of  the  fluid  past  the  rudder 

V  ship  velocity 

V.  maximum  ship  velocity 

W  weight  of  rudder 

w  width  of  rectangular  slot 

X  valve  flapper  displacement 

y  hydraulic  actuator  displacement 

Z  Laplace-domain  spool  valve  displacement 

z  valve  spool  displacement 

z'  effective  valve  spool  displacement  after  compensation  for 

spool  dead  zone  effects 

z''         effective  valve  spool  displacement  after  compensation  for 
spool  dead  zone  and  multiple  flow  gain  effects 

a  rudder  angle  of  attack  to  oncoming  stream 

3  time  scaling  factor  for  analog  program 

^-C„„        rudder  turning  rate  -  dependent  portion  of  C^„ 

ML.  rudder  turning  rate  -  dependent  portion  of  the  rudder 

twisting  moment 

e  eccentricity  of  piston  in  sleeve 

V  kinematic  viscosity 

M  absolute  viscosity 

gU 
Q  da 


6  rudder  angle 

6  desired  rudder  angle 


A  superscripted  dot  (')  denotes  the  rate  of  change  of  the  quantity 
with  respect  to  time. 

The  subscripts   and     denote  maximum  values  of  the  subscripted 

^  ^ ,        m     max  '^ 

variable. 


I.   THE  NATURE  OF  THE  PROBLEM 

A.   BACKGROUND  INFORMATION 

In  1971,  the  Naval  Ship  Research  and  Development  Center  (NSRDC) ,  in 
Annapolis,  Maryland,  identified  a  potential  stability  problem  in  an 
electrohydraulic  servosystem  which  was  being  developed.   The  control 
system  was  being  designed  to  control  the  movement  of  a  submarine's  rudder 
in  accordance  with  commands  from  the  helmsman.   The  NSRDC  developed  an 
analog  computer  simulation  (hereafter  referred  to  as  the  "baseline  simu- 
lation") which  incorporated  many  of  the  characteristics  of  an  actual 
hydraulic  control  system,  including  friction,  hydraulic  and  mechanical 
characteristics,  servovalve  behavior,  and  interdependence  of  elements 
necessary  to  create  a  simulation  of  the  actual  system.   The  simulation 
was  based  on  equations  developed  to  model  an  existing  system  used  in 
the  USS  PARGO  (SSN  650) .   Provisions  were  made  so  that  the  simulation 
could  be  adjusted  by  changing  values  of  the  various  components  of  the 
model,  in  order  to  conform  to  possible  changes  of  the  parameters  of 
the  system  under  development,  such  as  changes  in  dimensions,  materials, 
and  operating  conditions.   The  provision  for  adjustment  of  the  simulator 
parameters  also  allowed  simulation  of  other  existing  hydraulic  systems 
of  this  type.   All  parameters  which  were  assumed  to  be  of  significant 
impact  on  the  performance  of  the  system  were  incorporated  into  the 
simulation. 

When  the  simulator  parameters  were  arranged  to  model  the  performance 
of  the  system  under  development,  the  model  did  not  conform  to  expected 
system  behavior;  oscillations  in  pressures  and  control  surface  turning 


10 


rates  were  observed  which  had  not  been  observed  in  similar  systems  in 
use.   Further  investigations  failed  to  reveal  any  obvious  errors  in  the 
simulation,  or  deficiencies  in  the  parameters  which  were  believed  to  be 
of  significant  impact  on  the  performance  of  the  system. 

This  report  is  a  result  of  a  project  initiated  in  support  of  the 
simulation  attempt  at  the  NSRDC.   The  purpose  of  this  project  was  to 
determine  the  reasons  for  the  observed  behavior  of  the  hydraulic  system 
model,  and  to  develop  appropriate  modifications  and  refinements  of  the 
baseline  simulation,  in  order  to  more  closely  model  actual  system 
performance. 

B.   SYSTEM  UNDER  STUDY 

1 .  Introduction 

The  following  paragraphs  describe  the  general  characteristics  of 
the  idealized  system  that  the  NSRDC  used  in  the  development  of  its  simula- 
tion.  The  description  does  not  completely  describe  any  real  system. 
Rather,  it  is  a  model  which  was  designed  to  include  all  of  the  steering 
system  characteristics  that  were  expected  to  be  of  significant  impact  on 
the  performance  of  the  system. 

2.  General  System  Descriptions 

The  steering  system  (Figure  1)  included  the  servovalve,  connecting 
tubing,  actuating  ram,  mechanical  linkages,  rudder,  and  minor  components. 
The  operation  of  the  system  is  described  below.   (Numbers  in  parentheses 
refer  to  corresponding  labeled  parts  on  Figure  1.) 

A  desired  change  of  rudder  angle  is  initiated  by  moving  the  control 
wheel  (1)  by  hand,  which  generates  an  electrical  signal  having  a  voltage 
proportional  to  the  degree  of  movement.   This  voltage  representing  the 
desired  control  surface  angle  is  then  compared  with  a  voltage  proportional 

11 
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to  the  actual  position  of  the  hydraulic  ram  (6).   An  error  signal  is 
generated  which  approximates  the  difference  between  the  desired  and 
actual  rudder  positions  (2).   This  signal  actuates  an  electric  torque 
motor  (3)  in  the  servovalve,  which  moves  a  small  flapper  (4)  located  in 
the  narrow  space  between  two  nozzles,  causing  an  increase  in  the  pilot 
pressure  on  one  nozzle,  and  a  decrease  on  the  other.   This  difference 
in  pressure  causes  the  spool  valve  (5)  to  move,  thereby  permitting 
pressure  and  subsequent  fluid  flow  to  be  applied  to  the  hydraulic  ram 
(6),  resulting  in  a  correcting  movement  of  the  control  surface  (7)  through 
the  mechanical  linkage.   The  value  of  the  error  signal  is  modified  by 
movement  of  the  hydraulic  ram,  and  the  entire  process  is  continuously 
repeated  until  the  error  signal  is  reduced  to  zero.   As  the  error  signal 
decreases,  the  servovalve  closes,  control  surface  movement  ceases,  and 
equilibrium  is  established  at  the  new  position. 
a*   Component  Descriptions 

(1)  Servovalve.   The  servovalve  (Figure  2)  used  in  a  submarine 
steering  control  system  is  a  two-stage,  force  feedback  model,  incorporating 
a  torque  motor  actuated  by  an  electrical  error  signal;  a  flapper-type 
pilot  valve,  the  flapper  itself  also  serving  as  the  armature  of  the  torque 
motor;  and  a  spool  valve,  actuated  by  the  fluid  flow  caused  by  the  movement 
of  the  pilot  flapper  valve.   Movement  of  the  spool  results  in  the  direction 
of  fluid  flow  to  the  actuating  ram. 

(2)  Actuator.   The  actuator  is  a  linear,  two-way  hydraulic 
ram.   Initially,  it  was  assumed  that  the  actuator  was  of  compensated  con- 
struction; that  is,  it  was  modeled  as  if  it  had  been  designed  with  equal 
effective  areas  on  both  sides  of  the  piston.   0-rings  are  incorporated 
into  the  design  of  the  actuator,  eliminating  essentially  all  leakage 
from  one  side  of  the  piston  to  the  other. 

13 
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(3)   Mechanical  Linkage.   The  linkage  connecting  the  hydraulic 
ram  to  the  rudder  consists  of  a  long  shaft,  support  bearings,  a  packing 
at  the  penetration  of  the  pressure  hull,  a  crosshead  at  the  end  of  the 
shaft  opposite  the  ram,  and  a  rod  connecting  the  crosshead  to  the  rudder 
tiller.   The  crosshead  provides  a  supported,  pinned  joint  from  which  the 
connecting  rod  originates.   It  was  assumed  to  be  constrained  from  lateral 
movement  by  an  open-ended  sleeve  which  was  flooded  with  sea  water. 

(A)   Rudder.   The  rudder  is  composed  of  two  portions,  an  upper 
rudder  and  a  lower  rudder,  interconnected  to  the  tiller  by  a  single  shaft, 
and  located  on  the  top  and  bottom  of  the  ship's  hull,  respectively.   The 
rudder  is  symmetric  along  the  chord  length,  and  is  designed  for  operation 
within  the  range  of  left  35  degrees  rudder  to  right  35  degrees. 

(5)  Other  Components.   Other  components  include  piping;  an 
electric  transmitter  associated  with  the  control  wheel;  an  error  feedback 
servo  associated  with  the  actuator;  a  hydraulic  pressure  supply  which  is 
assumed  for  design  purposes  to  be  constant  in  pressure,  and  not  subject 
to  influence  by  steering  system  conditions;  and  a  petroleum-base  hydraulic 
fluid. 

C.   NSRDC  BASELINE  SIMULATION 

The  following  sections  describe  the  work  done  at  the  NSRDC  in  order 
to  simulate  the  submarine  steering  control  system.   The  governing  equations 
for  each  of  the  components  are  stated  as  they  were  derived  by  the  NSRDC, 
along  with  any  necessary  assumptions.   A  schematic  circuit  diagram  of 
the  analog  computer  arrangement  is  presented,  and  comments  from  the 
NSRDC  about  the  nature  of  the  problem  and  suspected  origins  of  the  oscil- 
lation are  stated. 
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1.   Governing  Equations  and  Assumptions 

The  governing  equations  that  were  provided  for  each  of  the 
components  are  listed  verbatim.  Only  the  governing  equations,  assump- 
tions, and  values  for  constants  were  communicated  to  the  writer;  the 
procedure  by  which  the  expressions  were  derived  is  unknown.   See 
Appendix  A  for  explanations  of  the  notation  used, 
a.   Servo  Valve 

The  servo  valve  equations  include  expressions  for  the 
valve  flapper  displacement,  spool  valve  displacement,  nonlinearity  com- 
pensations, and  fluid  flow  rates  through  the  valve. 

(1)   Valve  Flapper  Displacement.   The  valve  flapper 
displacement,  x,  is  expressed  as 


sat 


-   (K-(6  -6)-K_z) 
K_    1   o      J 


+.003 
-.003 


(1) 


where  K  is  the  amplifier  gain  coefficient;  K-  is  the  flapper  centering 

spring  rate,  9  is  the  desired  position  of  the  rudder;  6  is  the  actual 
o 

position  of  the  rudder,  K  is  the  wand  spring  rage,  and  z  is  the  present 
valve  spool  position.  No  explicit  assumptions  were  given  with  regard  to 
this  equation. 

(2)  Valve  Spool  Displacement.   The  valve  spool  displacement 
is  obtained  through  the  variation  with  time  of  the  valve  spool  velocity  z: 


0.5   z  >.5  and  SGN(z)  »  SGN(x) 


K.x  otherwise 
4 


(2) 


where  K,  is  the  pilot-stage  valve  flow  coefficient.   This  expression  is 
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best  explained  by  stating  that  the  spool  velocity  is  proportional  to  the 
flapper  valve  displacement,  until  the  spool  is  fully  displaced  against 
the  valve  body,  where  the  movement  must  of  course  cease.   The  velocity 
again  becomes  non-zero  only  if  the  direction  of  movement  is  reversed,  and 
the  valve  is  allowed  to  travel  away  from  the  hard  restraining  block. 

Nonlinear  effects  in  the  valve,  such  as  varying  flow  gain 
(i.e.,  flow  is  not  proportional  to  spool  displacement,  due  to  the  shape 
of  the  ports  used  in  the  spool  valve),  spool  friction,  deadband,  and  torque 
motor  hysteresis  are  included  in  the  next  equations. 

Dead  band  is  included  in  the  system  through  the  implementa- 
tion of  the  following  equations: 

z  -  db   where   z  >db,  7>0 

.     I 

7.     =  <C       z  +  db   where   |zl  >db,  z<0  (5) 

0       if   -db<z<db 

where  z'  is  the  valve  spool  displacement  after  adjustment  for  spool  dead 
zone.   These  expressions  give  the  effective  spool  displacement,  i.e., 
adjust  the  spool  displacement  to  the  value  it  would  have  for  the  given 
flow  if  no  dead  zone  were  constructed  in  the  valve. 

Multiple  flow  gain  is  accommodated  by  combining  z'  with 


z"  = 


5z'  +  .5z'l  z'l  .  (4) 

where  z"  is  the  valve  spool  displacement  after  spool  dead  zone,  and  after 
multiple  flow  gain. 

Compensation  for  other  nonlinearities  such  as  spool 
friction  and  torque  motor  hysteresis  may  have  been  accomplished  through 
adjustment  of  the  coefficients  in  the  equations  for  spool  dead  band,  and 
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multiple  flow  gain.   No  additional  expressions  for  these  nonlinearities 
were  given. 

(3)   Flow  Equation.   The  flow  rate  in  the  spool  is  related  to 
the  corrected  spool  displacement  z",  the  load  pressure  drop  P, ,  and  the 
supply  pressure  P  by  the  following  expression: 


s 


Q  =  z"K3Sgn(P3-sgn(z)PL)y  Pg  -  sgn(z)PL  (5) 

where  K^  is  the  main-stage  valve  flow  coefficient.   No  explicit  assumptions 
were  listed  for  this  equation, 
b.   Actuating  Ram 

The  position  of  the  actuating  ram  is  described  in  two  ways;  by 
the  position  of  the  rudder  itself,  and  by  the  integral  of  the  flow  rate 
which  has  passed  through  the  valve. 

(1)  Ram  Position  in  Terms  of  Rudder  Position.   The  ram  position 

y  is  related  to  the  rudder  position,  with  correction  for  mechanical  resilience 
in  the  linkage,  by 

APr 

y  =  Lsine   +  -jt —  .  (6) 

where  L  is  the  length  of  the  tiller,  9  is  the  rudder  angle  in  degrees,  K 

is  the  mechanical  resilience  coefficient  of  the  system,  and  A  is  the  actuator 
piston  area. 

(2)  Ram  Velocity  Expressed  in  Terms  of  Pressure  Change  and  Flow. 
The  NSRDC  determined  that  the  ram  velocity,  flow  rate,  and  load  pressure 
have  the  following  relationship  to  one  another: 

Q  =  -Ay  4-  K^P^  (7) 

where  K  is  the  system  hydraulic  resiliance,  and  Pj  is  the  rate  of  change 
6  ^ 

of  load  pressure  drop. 
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In  the  development  of  these  equations  it  was  assumed  that 
there  was  no  leakage  from  the  control  valve  to  the  actuating  cylinder  at 
line  connections,  seals,  and  so  forth,  and  that  there  was  no  inLcrnal 
leakage  in  the  actuating  cylinder  itself.   Also,  it  was  assumed  that 
the  actuator  used  was  compensated,  or  incorporated  equal  areas  on  opposite 
sides  of  the  piston. 

c.   Forces  and  Torques  on  the  System 

All  forces  on  the  system  were  expressed  in  terms  of  the  equiva- 
lent torques  which  were  placed  on  the  rudder  structure,  and  a  moment 
balance  was  written: 

le  =  T  +  T„  +  T_  (8) 

p     H     F 

where  1  is  the  polar  moment  of  inertia  of  the  rudder  mass  and  the  included 

water  mass,  6  is  the  angular  acceleration  of  the  rudder,  T  is  the  torque 

P 

developed  as  a  result  of  the  pressure  drop  across  the  actuator,  T  is  the 

hydrodynamic  torque,  and  T  is  the  friction  torque. 

r 

(1)   Friction  Torque.   T  is  the  most  complicated  of  the 
torques  to  describe.   It  is  the  combined  friction  torques  from  the  shaft 
seals,  the  resultant  of  the  normal  force  on  the  rudder  bearing  surface, 
and  any  other  sources  of  sliding  (coulomb)  friction  present  in  the  system. 
Stationary  objects  in  contact  with  each  other  stick  together  until  sub- 
jected to  forces  larger  than  the  friction  force  present  between  the 
objects.   Additionally,  the  level  of  applied  force  necessary  to  sustain 
movement  at  very  low  speeds  usually  exceeds  the  force  needed  at  higher 
velocities.   In  order  to  model  these  system  characteristics,  the  NSRDC 

developed  the  following  expi^assions  for  T  : 

r 
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Tp   = 


-Tssgn(VTp)    if   |VT^>Tsj 
-T^sgn(e)  lel>e 


(9) 


where  e  is  theoretically  zero,  but  in  practice  is  set  to  approximately 
0.2  percent  of  the  maximum  angular  velocity  of  the  rudder.   T  =  kTg , 
where  k  is  slightly  less  than  one,  typically  0.95,  and  T„  is  expressed  by 


Tg  -  Cj^R 


W  +  /(Kge)2+Kjje+A-PL) 


+  K 


10 


(10) 


where  C^  is  the  coefficient  of  friction  of  the  stock  bearing,  W  is  the 
weight  of  the  rudder,  R  is  the  rudder  stock  radius.  Kg  is  the  lift  co- 
efficient of  the  rudder,  K-, -,  is  the  drag  coefficient  of  the  rudder,  and 
K,|^  is  the  seal  friction  coefficient.   It  should  be  noted  that  the 

expression  for  T  does  not  include  any  terms  for  damping  proportional 
f 

to  rudder  turning  rate, 

(2)   Hydrodynamic  Torque.   T  is  the  torque  generated  by  the 
rudder  when  it  is  placed  at  some  angle  of  attack  other  than  zero,  due  to 
the  fact  that  the  center  of  pressure  on  a  rudder  does  not  coincide  with 
the  turning  axis  of  the  structure.   The  hydrodynamic  torque  is  expressed 


as 


^H=(y^S^-V^'  ) 


(11) 


where  V  is  the  ship  velocity,  V  is  the  maximum  ship  velocity,  and  Ky  and 

K  are  hydrodynamic  torque  load  coefficients.   By  adjusting  the  values  of 
8 

Ky  and  Kg,  the  hydrodynamic  torque  variation  with  rudder  position  can  be 
adjusted  to  approximate  almost  any  rudder  design.   If  K  and  Kg  are 
redefined,  the  above  expression  can  be  written  in  nondimensional  form: 
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FIGURE    3 
HYDRODYNAMIC    TORQUE    BEHAVIOR 
FOR    VARIOUS    VALUES   OF  Ky    AND    Kq 
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"^Hmax    '^  \  max/    ^  l  max 


max 


(12) 


where  T^^j^^^  is  the  maximum  hydrodynamic  torque  corresponding  to  Vq,  and 
®max  ^^  ^^^^   maximum  rudder  angle.   See  Figure  3  for  curves  of  Tj,/T 
for  varying  values  of  K,  and  K  .   It  can  be  seen  that  a  variety  of  rudder 
torque  characteristics  can  be  simulated  by  adjustment  of  K^  and  K„. 

I  o 

(3)  Pressure  Torque.  The  pressure  torque  Tp  is  the  torque 
that  is  placed  on  the  stock  of  the  rudder  by  the  hydraulic  ram,  through 
the  connecting  linkage.   It  can  be  expressed  as 

T  =  LAP^cos  e  ,  (13) 

where  L  is  the  tiller  length,  A  is  the  area  of  the  actuator  piston,  and 
Pj  is  the  load  pressure  drop. 

2.  Analog  Computer  Simulation 
a.   Hardware  and  Parameters 

Figure  4  is  a  schematic  diagram  of  the  analog  computer  program 
which  was  designed  to  model  the  governing  equations.   Table  I  gives  a  list 
of  typical  system  dimensions  and  other  parameters  which  were  used  to  set 
simulation  values.   Table  II  gives  a  list  of  the  potentiometer  settings 
for  the  values  listed  in  Table  I,  and  also  lists  the  expressions  which 
were  used  to  develop  the  settings. 

3.  Discussion  of  the  Baseline  Simulation 
a.   Performance  Curves 

Examination  of  sample  recordings  obtained  from  the  NSRDC 
showed  that  oscillations  in  velocity  and  pressure  were  present  in  the 
simulation  outputs.   The  NSRDC  indicated  that  the  oscillation  appeared 
to  be  a  function  of  the  inertia  of  the  steering  linkage,  ram,  and  rudder, 
and  of  the  system  mechanical  and  hydraulic  resiliences.   Variation  of 
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Table  I 
Submarine  Steering  System  Parameters  Used  in  the  Baseline  Simulation 

Parameter  Value  Parameter  Value 

A  80  in^  L  30  in 

Cd  -05   ■  Pmax  6,000  psi 

db  .05  P,  3,000  psi 

Lmax  '    ^ 

I  9100  in-lbf-sec^    Pg  3,000  psi 

K  .6   .  Qj^^  278,86  in-^/sec 

K2  110  R  18  in 

K^  5.0  Tpjj,^  l.ieeiXlO^in-lbf 

K.  375  T„  3.93X106in-lbf 

4  Hmax 

Kr  7.20  T^  7.2X106in-lbf 
-^                          Pmax 

K^  0.015  W  25,000  Ibm 
0 

K_  1.5  X  .002  in. 
'                          max 

Ko  2.5  y  17  in. 

8  -^max 

K„  30,000  z  .5  in. 

9  max 

KjQ  150,000  ^^^^  35  degrees 

K^^  3,000  \^  10  deg/sec. 

K,2  800,000  6,  time  scaling  factor  1 
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Table  II 

Potentiometer  Expressions  and  Settings  for  Baseline  Simulation 

(See  Fig.  4) 


Potentiometer 
No. 


6 
7 
8 
9 
10 

11 

12 

13 
14 

15 


Expression 


lOOK^x 

2  max 


Simulation 
Technique 

'^4^max 


lOz 


max3 


%^max 

lOOK^x 

2  max 


db 
db 

multiple  flow  gain 
multiple  flow  gain 
^max     (1st  footnote) 

WY^ 

AL%ax(TT/180)C2nd    footnote) 
20  K^  PLmax  B 


^Lmax/p 

^^/Pmax 

^^Lmax 
I2-'max 


max 


Setting 
.9546 

.9546 

.1 
.1500 

.1136 

.05 

.05 

.5 

.5 

.406 

.4654 

.5000 


0177 


8824 


2y 


max 
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Potentiometer 

No. 


16 

17 
18 
19 

20 

21 
22 


Expression 

AP,        (3rd  footnote) 
3.45X105 

^ll^max    ^^^^  footnote) 
3.45X105 

'^Q^max     ^■^^'^    footnote) 
1.075X106 


Setting 


6957 


3044 


9768 


W       (3rd  footnote)  Q233 

1.075x106 

7  AcvTn^  ^(3rd  footnote) 

r3. 45X10^/        ^.2  .6934 

^        /l. 129x106)^  •"''■^^ 

r  ^ (3d  footnote) 

(1.075X106/^_^29X10^)2  .9066 


KlO 
Fmax 


.1286 


23 


Cd  R(1.129X106)     (3rd  footnote)   .8714 

T 
Fmax 


24 
25 


Tp 


max 


looe^^axi 


.0010 
.1282 
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TFmpY(Tc/Ts) 

lOOe    13 
max 


.1269 


27 


max 


'max 


2857 


28 


29 


max     (4th  footnote) 


200 


10 


1750 


25 


26 


Potentiometer 

No. Expression  Setting 


30  K7_ 

10 


15 


31  '^Hmax  .3370 


32  Tpmax  .6175 

10  Tpniax 


1.  This  expression  was  later  found  to  be  incorrect.   See  the  text, 
Initial  Investigation. 

2.  This  expression  was  later  found  to  be  incorrect.   See  the  text. 
Initial  Investigation. 

3.  Numbers  appearing  in  the  numerators  or  denominators  of  these  expressions 
represent  scaling  values  calculated  by  the  NSRDC. 

4.  Adjustment  for  SIN  -  COS  generator  design  requirements. 
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these  parameters  modified  the  fundamental  frequency  of  the  oscillation, 
but  did  not  eliminate  it.   It  was  found  that  the  oscillation  could  be 
eliminated  by  the  artificial  introduction  of  viscous  damping  in  the 
system,  but  only  in  much  greater  amounts  than  seemed  to  be  present  in 
the  actual  system. 

b.   Suspected  Origins  of  the  Oscillations 

After  some  study  of  the  baseline  simulation,  the  NSRDC  folt 
that  the  problem  was  located  in  either  their  recording  equipment  associated 
with  their  analog  computer,  in  the  mathematical  models  which  had  been 
developed,  in  settings  of  the  simulator  parameters,  or  in  the  basic  assump- 
tions used  in  the  derivations  of  the  governing  equations. 
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II.   INITIAL  INVESTIGATION 

The  initial  investigation  was  conducted  in  order  to  gain  a  more 
complete  understanding  of  the  simulation,  and  to  determine  problem  areas 
requiring  extensive  investigation.   The  baseline  simulation  was  duplicated 
on  a  large  analog  computer,  techniques  and  models  currently  used  in  simulat- 
ing hydraulic  control  systems  were  studied,  and  information  was  gathered 
from  personnel  familiar  with  the  baseline  simulation  development  and  sub- 
marine steering  control  systems.   The  baseline  simulation  was  continually 
evaluated  in  terms  of  the  information  which  had  been  gained,  and  several 
problem  areas  were  found  which  required  further  analysis. 

A.   DUPLICATION  OF  THE  BASELINE  SIMULATION 

The  baseline  simulation  required  duplication  on  local  equipment  in 
order  to  provide  a  check  of  the  correctness  of  the  programming  completed 
at  the  NSRDC.   The  Naval  Postgraduate  School's  COMCOR  CI -5000  analog 
computer  was  used  in  this  effort.   The  CI -5000  is  a  general  purpose 
computing  system,  including  analog,  digital  control,  and  patchable  logic 
sections,  and  incorporates  enough  components  of  various  types  to  accommo- 
date the  baseline  simulation. 

The  writer  first  attempted  to  directly  adapt  the  baseline  simulation 
analog  schematic  directly  to  the  CI-5000,  adjusting  the  simulation  only 
to  compensate  for  CI-5000  peculiarities,  thereby  assuming  that  the  equations 
and  potentiometer  expressions  of  the  baseline  simulation  were  correct,  and 
that  the  analog  schematic  diagram  accurately  modeled  the  governing  equations 
This  first  attempt  was  unsuccessful;  analog  components  rapidly  overloaded, 
and  the  system  failed  to  reproduce  the  sample  curves  that  the  baseline 
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simulation  at  the  NSRDC  had  produced  when  operated  under  identical 
conditions.   Consequently,  the  baseline  simulation's  governing  equations, 
potentiometer  expressions,  and  analog  programming  were  checked  for  mathe- 
matical and  functional  correctness. 

The  governing  equations  were  found  to  be  consistent  with  the  sign 
conventions  evidently  adopted  by  the  developer  (See  Figure  1),  with  the 
exception  of  equation  (7), 

Q  =  -Ay  +  K.P,  . 

Examination  of  the  geometry  of  the  system  (Figure  1)  revealed  that  the 
sign  of  the  y  -  term  should  have  been  positive.   Thus,  equation  (7)  should 
have  read 

Q  =  Ay  -(-  K.P,  .  (7a) 

O  L 

The  NSRDC  later  confirmed  that  this  was  the  case. 

Potentiometer  expressions  were  found  to  be  consistent  with  the 
baseline  simulation  schematic  diagram  and  the  governing  equations,  witii 
the  exception  of  potentiometers  10  and  11,  which  were  associated  with 
integrator  9  of  the  baseline  simulation  (Figure  4).   A  term  in  the  denom- 
inator of  each  potentiometer  expression  apparently  had  been  dropped;  the 

2 
denominators  should  have  read  10 (P,    )[K^-«-(A  /K, .)]  instead  of  10 (KP,    ) 

Lmax   6      12  6  Lmax 

and  20(K-P,    )  for  pots  10  and  11,  respectively.   The  need  for  the 
o  Lmax 

additional  term  in  each  expression's  denominator  was  later  confirmed  by 
the  NSRDC.   The  additional  terms  modified  the  values  of  the  two  potenti- 
ometers significantly;   the  value  of  pot  10  changed  from  0.6197  to  0.4060; 
the  value  of  pot  11  changed  from  0.4654  with  a  gain  of  20  on  integrator 
9  to  0,6060  with  a  gain  of  10.   Tests  later  showed  that  the  omission  of 
these  coefficient  terms  was  the  primary  reason  for  the  initial  failure 
to  duplicate  the  baseline  simulation. 
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Except  for  the  two  potentiometer  settings,  the  analog  program  (Fig.  4) 
appeared  to  accurately  represent  the  verified  governing  equations.   Accord- 
ingly, another  attempt  was  made  to  adapt  the  program  to  the  CI-5000.   Only 
slight  deviations  from  the  baseline  schematic  diagram  were  necessary  in 
order  to  accommodate  the  CI-5000  component  peculiarities. 

When  the  adapted  baseline  simulation  was  operated  under  the  potenti- 
ometer conditions  listed  in  Table  II,  with  the  exception  of  pots  9  and  10 
which  were  set  as  previously  modified,  the  system  behaved  in  the  same 
manner  as  it  had  during  operation  at  the  NSRDC.   Oscillations  identical 
to  baseline  forms  were  observed  in  6  and  P   (Figure  5).   As  a  result  of 
this  duplication  of  performance,  using  several  different  types  of  recording 
equipment  that  exhibited  varying  response  characteristics,  it  appeared  that 
the  recording  equipment  at  the  NSRDC  was  not  the  source  of  the  oscillations. 
Further,  since  all  changes  in  mathematics  and  in  the  settings  of  the  simu- 
lator parameters  were  implemented  before  the  baseline  performance  was 
successfully  duplicated,  it  appeared  that  the  errors  which  were  found  were 
generated  through  the  communication  of  the  baseline  simulation  to  the 
writer,  and  were  not  present  in  the  program  used  to  generate  the  sample 
curves.   Apparently,  then,  the  mathematical  models  were  consistent  in 
sign  conventions  and  had  been  manipulated  correctly  in  the  development  of 
the  analog  board,  the  potentiometers  were  correctly  scaled,  and  the  analog 
wiring  was  correct.   Thus,  the  simulation  appeared  to  be  correct  within 
the  framework  of  the  fundamental  assumptions  involved  in  the  mathematical 
formulations;  however,  the  oscillations  were  still  present  in  the  output. 

Two  situations  seemed  to  be  possible.   First,  perhaps  the  oscillations 
evidenced  in  the  NSRDC  data  were  not  representative  of  the  behavior  present 
in  real  submarine  steering  systems,  and  the  simulation  was  therefore 
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inadequate.   This  situation  could  have  occurred  because  omissions  were 
made  during  the  formulation  of  the  baseline  simulation  governing  equations, 
such  as  a  miscalculation  of  the  significance  of  some  physical  process 
present  in  a  steering  system.   The  situation  could  also  have  occurred 
because  system  materials,  dimensions,  and  properties  were  incorrectly 
represented  in  the  development  of  the  numerical  values  using  in  the  "typical" 
system  simulation,  due  to  the  inclusion  of  inconsistent  units  conversions, 
unrealistic  dimensions,  or  faulty  physical  materials  constants. 

Second,  perhaps  the  oscillations  appearing  in  the  NSRDC  data  were 
representative  of  real  system  behavior,  and  the  difficulty  had  its  origins 
in  the  harboring  of  an  incorrect  knowledge  of  the  actual  performance  of  a 
ship's  steering  control  system.   No  actual  system  recordings  were  submitted 
with  the  baseline  simulation,  and  therefore  no  immediate  comparison  was 
possible  between  the  simulation  performance  and  real  conditions  in  the 
system. 

B.   BACKGROUND  INVESTIGATION 

1.   Study  of  Currently  Used  Models 

Assuming  that  the  oscillations  observed  in  the  basdline  simulation 
were  not  representative  of  real  steering  system  behavior  led  to  a  back- 
ground investigation  of  techniques  and  models  currently  used  in  simulating 
hydraulic  systems  in  order  to  develop  a  clear  knowledge  of  the  fundamental 
assumptions  made  during  the  development  of  the  baseline  simulation.   This 
background  examination  consisted  of  studies  in  the  hydraulic  controls 
field,  and  communication  with  personnel  familiar  with  the  simulation  and 
ship  steering  control  systems. 

The  search  for  currently  used  models  and  techniques  produced  several 
results.   It  was  noted  that  if  the  expressions  for  the  baseline  simulation 
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were  transformed  into  Laplace  notation,  and  deadband  and  multiple  flow 
gain  effects  were  temporarily  neglected,  the  resultant  expression  would 

be 

K  /K 
Z(s)  =   ~±~^ —  E(s)  (lA) 

s-hK3/K2 

where  Z  and  E  represent  the  Laplace-domain  variables  for  z,  the  valve 
spool  position,  and  e,  the  error  signal,  respectively,  and  s  is  the 
Laplace  transform  variable.   This  function  was  identical  in  form  to  one 
developed  by  Morse  ( 1  J  for  a  valve  incorporating  centering  springs  instead 
of  force  feedback,  for  use  in  systems  where  the  load  was  assumed  Lo  hi. 
very  small  compared  with  the  available  hydraulic  power.   The  actual  valve 
being  used,  however,  consisted  of  a  two-stage,  force-feedback  design  being 
used  to  control  a  load  of  the  same  order  of  magnitude  as  the  available 
hydraulic  power.   Merritt  [2]  developed  a  transfer  function  block  diagram 
for  a  two-stage,  force-feedback  servovalve,  which  incorporated  the  general 
form  shown  in  Figure  6.   A  vast  degree  of  simplification  of  Merritt's 
expression  obviously  would  have  been  necessary  in  order  to  arrive  at  the 
expression  used  in  the  baseline  simulation.   Quite  possibly,  the  coef- 
ficients indicated  by  Merritt  in  many  instances  were  of  such  magnitudes 
that  their  effects  were  not  important  in  the  actual  performance  of  sub- 
marine steering  system  servovalves,  since  high  frequency  effects  would 
become  unimportant  when  other  system  frequencies  were  quite  low.   Further- 
more, the  pressure  feedback  to  the  flapper  valve  that  Merritt  included  in 
his  expression  was  possibly  unimportant  in  this  particular  case.   However, 
it  seemed  advisable  to  check  the  numerical  value  of  each  term  in  the 
detailed  expression,  in  order  to  evaluate  the  validity  of  the  baseline 
simulation's  expression. 
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It  was  found  [3]  that  many  of  the  servovalve  design  values,  needed  to 
calculate  the  coefficients  in  the  Merritt  block  diagram,  were  difficult 
r.o  obtain.   It  was  also  found  that  the  expression  used  in  the  baseline 
simulation  was  considered  to  be  sufficiently  accurate  for  the  particular 
design  studies  in  which  it  would  be  used.   Finally,  it  was  determined 
through  observations  of  the  performance  of  the  adapted  simulation  program 
that  the  oscillations  did  not  appear  to  be  affected  by  the  servovalve 
action.   Thus,  further  investigation  of  the  validity  of  the  servovalve 
simulation. did  not  appear  to  be  justified  until  the  origin  of  the  oscilla- 
tions had  been  fully  investigated. 


SERVO 
VALVE 


ACTUATING 
RAM 


RUDDER  & 
LINKAGE 


Figure  7.   Functional  Block  Diagratr.  of  a  Submarine  Steering  System 

Figure  7  is  a  function  block  diagram  of  a  submarine  steering 
control  system.   The  control  system  was  not  entirely  closed  loop  in  con- 
struction, since  the  error  «%ignrii  controlling  the  function  of  the  servo- 
valve was  generated  by  the  ram  position  rather  than  by  the  actual  rudder 
position.   This  placement  of  the  feedback  element  was  consistent  with 
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established  design  practice  [4],  since  the  inclusion  of  the  compliant 
structure  of  the  rudder  and  linkage  within  the  servo  loop  could  result  in 
significant  instability  problems.   Additionally,  an  electrical  feedback 
element  located  outside  the  pressure  hull  of  a  submarine,  in  a  location 
flooded  with  sea  water,  would  pose  several  design  problems,  including  the 
need  for  increased  reliability  due  to  its  inaccessability  when  at  sea. 
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Figure  10.   Circuit  Modification  of  Fig.  4, 
Introduction  of  Viscous  Damping 


Examination  of  the  curves 
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(Figure  8)  generated  by  the  CI-5000  simulation  revealed  that  the  oscilla- 
tions were  present  only  in  the  linkage  and  the  load  portions  of  the  simu- 
lator.  Manipulations  of  the  baseline  simulation  friction  coefficients 
failed  to  eliminate  or  modify  the  effect,  although  total  exclusion  of  any 
friction  in  the  system  resulted  in  only  very  small  undamped  oscillations 
in  9,  and  noticeable  oscillations  in  6  and  P  after  the  control  surface 
had  essentially  reached  the  desired  position  (see  Figure  9).   Further 
tests  showed  that  small  values  of  negative  feedback  across  integrator  25, 
which  were  introduced  through  the  use  of  the  circuit  modification  indicated 
in  Figure  10,  resulted  in  elimination  of  the  oscillations.   This  feedback 
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across  the  integrator  would  be  representative  of  damping  in  an  actual 

system  containing  such  elements  as  dashpots.   With  no  other  friction  in 

the  system,  a  feedback  level  of  0.005  -r—  was  sufficient  to  eliminate  the 

limit  cycle.   Furthermore,  feedback  levels  of  0.12  t—  were  sufficient  to 

6m 
approach  a  critically  damped  condition  in  the  load  and  linkage  simulator. 

(See  Figures  11  and  12.)   If  friction  levels  present  in  the  baseline 

simulation  were  programmed  into  the  CI-5000,  approximately  the  same  value 

of  B  -  damping  as  in  the  frictionless  configuration  was  required  to  produce 

critically  damped  behavior,  while  the  sustained  limit  cycle  was  eventually 

eliminated  by  the  presence  of  the  baseline  friction  itself. 

Damping  terms  proportional  to  v   were  totally  absent  from  the  baseline 
simulation  equations.   Intuitively,  however,  some  damping  would  be  present 
in  any  system  in  which  viscous  fluids  such  as  sea  water  are  displaced  by 
moving  objects,  since  work  must  be  expended  to  keep  objects  moving  at  a 
constant  rate,  whether  linearly  or  rotationally,  through  a  real  fluid. 
Thus,  the  neglect  of  damping  in  the  linkage  and  rudder  portion  of  the 
simulation  did  not  appear  to  be  realistic.   A  superficial  examination  of 
the  system  did  not  reveal  any  large  sources  of  damping,  or  even  sources 
which  would  provide  a  fraction  of  that  needed  to  approach  a  critically 
damped  condition  in  the  linkage.   Anticipating  that  an  investigation  of 
system  damping  could  become  quite  involved,  the  writer  identified  this 
area  as  requiring  first  priority  for  further  investigation,  and  postponed 
further  analysis  of  it. 

Examination  of  steering  control  system  drawings  for  the  SSN-671  class 
submarine  revealed  that  functional  concepts  of  several  features  of  the 
design  had  been  altered  in  the  process  of  the  development  of  the  idealized 
model.   A  compensated  actuator  had  been  assumed  in  the  simulation  devel- 
opment; i.e.,  equal  areas  had  been  assumed  on  both  sides  of  the  actuator 
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ram  (Figure  1),  while  examination  of  the  drawings  revealed  that  the 
actuator  on  a  typical  submarine  was  usually  uncompensated.   It  was 
learned  that  the  effect  of  neglect  of  this  departure  from  actual  system 
design  was  thought  to  be  small,  and  not  important  in  terms  of  the  projected 
simulation  use  [5). 

Further  study  of  the  drawings  revealed  that  the  crosshead  configu- 
ration had  been  simplified.   The  crosshead,  originally  assumed  to  operate 
in  a  sleeve  open  at  both  ends.  Figure  1,  was  in  reality  supported  by  a 
sleeve  open  at  only  one  end  in  most  submarine  designs.  Figure  13.   Further 
investigation  of  the  crosshead  influence  on  the  system  seemed  to  be  justi- 
fied, since  the  assembly  resembled  a  dashpot  in  appearance,  and  possibly 
was  capable  of  providing  the  source  of  damping  which  seemed  to  be  needed 
in  the  system.   Consequently,  a  detailed  analysis  of  the  effect  of  a 
typical  crosshead  structure  on  steering  system  performance  seemed  to  be 
in  order,  and  was  identified  as  an  area  for  further  work. 
2.   Verification  of  Typical  System  Values 

Several  nondimensional  coefficients,  load  parameters,  and  properties 
had  been  used  to  develop  the  simulation,  such  as  the  hydraulic  resilience 
coefficient,  the  rudder  mass,  and  inherent  materials  properties.   Inaccurate 
development  of  system  numerical  values  could  have  contributed  to  the 
system  oscillations.   Attempts  were  made  to  determine  the  values  of  many 
of  these  constants  from  basic  principles,  but  the  effort  was  thwarted 
due  to  a  lack  of  information  regarding  dimensions,  physical  constants, 
system  construction  details,  and  algebraic  expressions  for  any  of  the  non- 
dimensional  constants  K  through  K^.   This  possible  source  of  inaccurate 
system  simulation  was  identified  for  further  investigation,  and  was  set 
aside. 
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3.   Actual  System  Performance 

Although  actual  system  recordings  of  e»  e,  and  P  were  mentioned 
16)  as  not  demonstrating  the  oscillations  found  in  the  simulation,  such 
recordings  were  not  forwarded  to  the  writer.   Because  an  accurate  simula- 
tion of  a  system  could  not  be  claimed  until  it  had  been  favorably  compared 
with  the  actual  system  performance,  it  was  felt  that  such  a  comparison 
would  be  beneficial,  and  the  subject  was  identified  as  requiring  further 
investigation. 

C.   SUMMARY  OF  THE  INITIAL  INVESTIGATION 
1.   Conclusions 

As  a  result  of  the  duplication  of  the  baseline  simulation,  and 
subsequent  experiments  performed  with  the  program,  the  following  conclu- 
sions were  reached: 

(1)  The  recording  equipment  at  the  NSRDC  had  accurately  recorded 
the  performance  of  the  simulation,  and  was  not  the  source  of  the  oscil- 
lations. 

(2)  The  equations  and  potentiometer  settings  were  consistent 
within  the  framework  of  the  fundamental  assumptions  involved  In  the 
simulation's  mathematical  development. 

(3)  The  oscillations  were  confined  to  the  load  portion  of  the 

simulator,  and  did  not  appear  to  be  sustained  by  overall  system  instability. 

L 

(4)  A  negative  feedback  level  of  0.005  -j —  was  sufficient  to 

.^m 

eliminate  the  limit  cycle  oscillation,  and  0.12 —  achieved  a  critically 

e 

damped  condition,  in  the  load  simulation  under  otherwise  frictionless 

n 

conditions.   A  negative  feedback  level  of  0.12  -; — was  also  sufficient  to 

e 

m 
create  a  critically  damped  condition  in  the  load  simulator  under  baseline 

friction  conditions. 
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Studies,  and  conununications  with  personnel  familiar  with  the 
simulation  and  with  actual  submarine  steering  control  systems,  led  to 
the  following  conclusions: 

(1)  The  servovalve  equation  used  in  the  baseline  simulation  was 
not  designed  to  completely  model  either  the  type  of  servovalve  used  in 
the  real  system  under  study,  or  to  be  used  with  the  type  of  load  present 
in  the  actual  system.   The  net  effect  of  the  substitution  of  a  simple 
expression  for  a  complicated  one  was  not  immediately  determined,  except 
that  it  did  not  appear  to  be  capable  of  affecting  the  observed  simulator 
load  oscillations. 

(2)  The  actuator  had  been  modeled  in  the  system  simulation  as 

a  compensated  ram,  rather  than  being  modeled  as  the  uncompensated  actuator 
used  in  real  systems.   Although  the  effect  of  this  simplification  on  the 
system  was  unknown,  communication  with  experienced  personnel  seemed  to 
indicate  that  the  effect  was  unimportant. 

(3)  The  actual  crosshead  design  had  been  modeled  in  such  a  way 
that  damping  effects  possibly  present  in  the  component  were  not  taken  into 
account.   Further  study  seemed  warranted. 

(4)  Viscous  damping  effects  in  the  load  had  been  neglected  in 
the  development  of  the  simulation,  although  such  effects  controlled  the 
unwanted  system  oscillations,  and  normally  would  be  present  to  some  degree 
in  any  system  incorporating  the  movement  of  structures  through  real  fluids. 
A  search  for  viscous  damping  sources  in  the  load  was  assigned  first  priority 
for  further  investigation. 

(5)  The  origins  of  many  of  the  nondimenslonal  constants  used  in 
the  simulation  were  not  obvious,  and  efforts  to  verify  their  correctness 
were  hindered  by  the  lack  of  their  related  algebraic  expressions. 
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(6)   The  realism  of  the  simulator  performance  needed  to  be 
checked  through  the  comparison  of  simulator  and  actual  system  recordings, 
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in.   LOAD  AND  LINKAGE  DAMPING 

Of  the  areas  identified  in  the  initial  investigation  as  requiring 
further  research,  only  the  search  for  sources  of  linkage  and  load  damping 
was  pursued  in  the  present  investigation.   Analysis  of  the  submarine 
system  identified  three  locations  where  damping  might  possibly  be  present: 
the  system's  dynamic  seals  and  bearings,  the  crosshead  structure,  or  the 
rudder  itself.   No  other  components  appeared  to  be  designed  in  such  a 
manner  that  they  could  possibly  produce  an  appreciable  quantity  of  damping. 

In  order  to  determine  levels  of  damping  in  the  various  components, 
detailed  information  was  needed  concerning  dimensions,  clearances,  materials, 
constructions,  and  lubrication  for  the  linkage  parts  and  rudder.   The  writer 
was  referred  [7]  to  Mr.  Leeland  Smith  [8]  of  the  Mare  Island  Naval  Shipyard, 
who  provided  detailed  drawings  of  the  steering  control  system  for  the  SSN 
671  class  submarine. 

A.   BEARINGS  AND  DYNAMIC  SEALS 

The  possibility  was  considered  that  the  expressions  employed  for  the 
friction  present  in  the  bearings  and  dynamic  seals  (equations  9  and  10) 
were  not  representative  of  the  actual  processes  present  in  the  submarine 
system,  and  that  the  inherent  friction  of  the  bearings  and  seals  was 
better  described  as  displaying  viscous  behavior. 

Keller  [9]  indicated  that  compression  packings,  such  as  those  used 
at  the  pressure  hull  penetration  and  the  actuator-piston  rod  interface, 
are  normally  lubricated  by  allowing  very  slight  amounts  of  leakage  past  the 
seal.   Examination  of  the  SSN-671  drawings  revealed  that  the  pressure  hull 
packing  was  additionally  lubricated  with  grease.   Despite  these  indications 
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that  viscous  effects  could  be  present  in  these  seals,  discussions  with 
Mr.  Smith,  William  Kendall  [lo]  ,  and  Leonard  Chase  [n]  of  the  Mare  Island 
Naval  Shipyard  resulted  in  the  writer's  learning  that  friction  in  such 
seals  was  best  approximated  as  sliding  friction,  since  friction  loads  were 
normally  observed  to  be  fairly  independent  of  actuator  rod  velocity,  due 
to  the  velocities  and  allowable  leakage  rates  present  in  the  system.   Thus, 
the  baseline  model  of  seal  friction  appeared  to  be  justified. 

Study  of  the  rudder  bearings  seemed  to  support  the  use  of  the  baseline 
simulation's  bearing  friction  expressions.   The  large  rudder  weight,  com- 
bined with  hydrodynamic  loads,  small  rates  of  rotation  present  in  the 
rudder  stock,  and  the  ill-defined  lubricating  qualities  of  sea  water, 
resulted  in  enough  analytical  uncertainties  being  present  to  discourage 
any  attempt  to  model  such  bearings  as  either  slipper-type  or  journal-type 
surfaces.   Again,  the  baseline  simulation  expression  seemed  to  provide  a 
bearing  surface  description  that  probably  was  as  accurate  in  magnitude  and 
parameter  dependence  as  any  other  expression  that  could  be  analytically 
developed  using  known  conditions  and  existing  lubrication  theory.   Experi- 
mental and  analytical  work  of  a  scale  beyond  the  scope  of  this  project 
could  prove  useful  in  future  design  work  involving  bearings  and  seals  of 
the  types  present  in  the  submarine  system. 

B.   CROSSHEAD 

1.   General  Description 

Examination  of  the  SSN-671  drawings,  and  discussions  with  Mssrs. 
Smith,  Kendall,  and  Chase  concerning  general  construction  details  of 
crossheads,  revealed  that,  as  a  rule,  crossheads  somewhat  resembled  dash- 
pots  in  construction,  due  to  their  pistqn-in  sleeve  arrangements,  and  their 
free-flooded  nature.   Provisions  were  normally  made  in  the  construction  of 
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crossheads  to  eliminate  any  possible  damping  effect  by  the  drilling  of 
large  holes  through  the  piston  portion  of  the  structure.   Because  damping 
sources  were  being  sought,  the  writer  decided  to  analyze  the  crosshead  in 
order  to  determine  what  damping  effects,  if  any,  were  present. 

Figure  13  shows  a  crosshead  design  based  on  information  taken  from 
the  SSN-671  drawings.   The  crosshead  consists  of  a  partially  bored  out 
cylinder  19  inches  in  length,  with  a  diameter  of  11.986  inches,  and  four 
one-inch  I.D.,  eight-inch  in  length  holes  axially  bored  through  the  solid 
head  of  the  piston.   This  piston  rides  in  a  twelve-inch  I.D.  tube,  typically 
sealed  along  its  length  except  at  the  rear  end,  which  is  open  to  allow  the 
passage  and  free  rotation  of  the  link  connecting  the  actuator  shaft  ex- 
tension and  the  rudder  tiller.   The  piston  incorporates  several  lock- 
preventing  0.75  in.  X  0.125  in.  circumferential  grooves,  and  two  0.125  in. 
X  0.5  in.  axial  slots  along  the  outside  radius.   As  a  result  of  the  radial 
clearance  between  the  piston  and  the  sleeve,  the  slots,  and  the  holes  in 
the  crosshead,  the  entire  sleeve  is  free  flooded,  and  flow  of  sea  water 
past  the  crosshead  must  take  place  as  it  is  moved  in  and  out  of  the  sleeve. 

2.  Assumptions 

In  approximating  the  behavior  of  the  flow  past  the  crosshead  and 
the  viscous  drag  associated  with  that  flow,  it  was  assumed  that  the  fluid 
velocities  through  the  holes,  slots,  and  radial  clearance  of  the  crosshead 
were  considerably  greater  than  the  motion  of  the  piston  relative  to  the 
sleeve,  thereby  implying  the  existence  of  fixed  orifices  of  irregular 
shape.   Laminar  flow  was  also  assumed. 

3.  Equations 

a.   Eccentric  Piston  in  a  Sleeve 
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For  an  annulus  formed  by  a  solid  eccentric  cylinder  in  a  sleeve 
with  no  relative  movement  involved,  the  flow  rate  Q  can  be  expressed  as  [12] 
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where  r  is  the  radius  of  the  sleeve, y  is  the  absolute  viscosity  of  the 
fluid,  1  is  the  length  of  the  tube,  e  is  the  eccentricity  of  the  piston  with 
respect  to  the  center  of  the  sleeve,  and  (P]^-P2)  is  the  pressure  drop  across 
the  length  of  the  tube. 

b.   Laminar  Flow  Through  Rectangular  Passages 
For  the  slots  [13] , 
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where  w  is  the  width  of  the  slot,  h  is  the  depths  and  all  other  symbols 

retain  their  previous  meanings. 

c.   Flow  Through  Short-Tube  Orifices 

Merritt  []_3J  ,  combining  the  work  of  Langhaar  [14]  ,  Kreith  [15]  , 

and  Shapiro  [15]  ,  developed  an  expression  for  the  flow  rate  Q^  through  short 

tube  orifices : 
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where  C  is  the  loss  coefficient,  p  is  the  fluid  density,  D  is  the  hole 
diameter.  Re  is  the  Reynolds  Number  based  on  the  mean  velocity  of  the 
flow  in  the  tube  and  the  hole  diameter,  and  all  other  variables  retain 
their  previous  meanings.   If  equations  15,  16,  and  17  are  manipulated 
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briefly,    the  pressure  drop  P^-P^   can  be   expressed   as  a  function  of  the  fl 
rate  Q^ ,    the  kinematic  viscosity  v,   D,   and  p: 


ow 
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If  a  total  mass  balance  is  made  across  the  crosshead,  it  is  found  that 

AcK  y  =  N^Q^  +  N^Q^  +  Q^  (20) 

where  A^^^^  is  the  cross-sectional  area  of  the  crosshead,  Nq  is  the  number  of 
holes  in  the  crosshead,  and  Ng  is  the  number  of  slots. 

The  force  caused  by  the  pressure  difference  across  the  crosshead 
can  be  expressed  as  an  equivalent  torque  T   at  the  rudder  pivot,  by  multi- 
plying the  crosshead  force  F   =  ^^1(^1-^2^    ^^  ^^^   length  L  of  the  tiller 
arm: 

^CH  =  A(.^  L  (P^-P2)  (21) 

Figure  14  shows  curves  of  T^^p^  versus  y  generated  through  the  solution  of 

equations  (18)  -(20)  using  the  previously  listed  crosshead  dimensions,  p 

3  c:    9 

equal  to  64  Ibm/ft  ,  and  v  equal  to  1.41  X  10'-^  ff^/sec  (water  at  50  degrees 

F) ,  for  various  sizes  of  piston  holes,  assuming  four  holes  and  two  slots, 

and  an  eccentric  piston  in  the  sleeve.   To  generate  the  curves,  equations 

(18)  and  (19)  were  solved  for  ^-,-^2   ^^  ^   function  of  Q-^.   Then  Q  and  (L. 

were  calculated,  and  appropriate  multiples  of  Q^,  Qg,  and  Q  were  used  in 

accordance  with  equation  (20)  to  produce  y.   Equation  (21)  was  then  used  to 

determine  T^., 
Ln . 

For  any  crosshead  incorporating  four  holes  any  larger  than  0.1 
inches  in  diameter,  it  was  found  that  zero  pressure  differential  existed. 
As  the  diameter  of  the  holes  was  decreased,  with  the  dimensions  and  number 
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of  slots  and  the  parameters  for  the  annulus  remaining  constant,  a  pressure 
differential  began  to  exist  across  the  structure.   Finally,  as  the  holes 
were  eliminated,  (^-^-^2^    approached  levels  sufficient  to  create  significant 
values  of  T^^^.      For  the  values  of  the  zero-holes  curve  used  in  Figure  14, 
'^CH  could  be  expressed  as 


T^^   =  7.48X103y  (22) 


y 

or,  since  -; — nj 


m    max 


TcH  =  (7.48X103:^ax)  (^^)  (23) 

If  the  above  equations  were  scaled  for  introduction  into  the 
simulation,  potentiometer  33  would  be  set  according  to  the  expression 

7.48X10^ymax  ,^.. 

a   =  ^-^ (24) 

33     100  le  max 
For  I  =  9100,  and  Xj^ax  -   ^   inches  per  second, 

^33   =  .00412. 

Since  this  value  of  the  potentiometer  represented  the  maximum  available 

Tp„,  it  was  seen  that  the  actual  effect  of  the  crosshead  on  the  system 

damping  was  very  small.   In  fact,  using  one-inch  diameter  holes  would  make 

the  crosshead  damping  essentially  zero.   However,  since  a  feedback  level  of 

.005  -» —  was  previously  found  to  be  sufficient  to  eliminate  limit  cycles  in 
^m 

the  "frictionless"  system  simulation,  .00412  ^—  could  assist  other  small 
sources  in  the  damping  of  the  oscillations  in  the  baseline  simulation.   Al- 
though the  results  of  the  investigation  of  the  crosshead  damping  were  less 
spectacular  than  had  been  desired  originally,  since  the  6-  feedback  sufficient 
to  create  a  critically  damped  condition  in  the  load  had  not  been  justified, 
the  author  felt  that  the  identification  of  a  small  source  of  damping  had 
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been  achieved,  that  identification  possibly  proving  to  be  of  some  worth  in 
future  design  work  if  a  small  amount  of  damping  somehow  proved  necessary. 
Some  argument  could  have  been  made  that  further  modification 
could  create  more  damping.   However,  for  this  to  be  done,  the  system  would 
need  to  be  sealed  at  both  ends,  due  to  the  approach  of  P2  to  the  vapor 
pressure  of  water  in  the  contained  volume  when  the  crosshead  was  moving 
out  of  the  tube.   Of  course,  the  open  end  of  the  crosshead  sleeve  was  not 
capable  of  being  closed,  due  to  the  primary  function  of  the  crosshead, 
which  was  to  .provide  a  portion  of  a  slider-crank  mechanism.   The  previously 
developed  levels  of  7^-,,  thus  appeared  to  be  the  maximum  practicable. 
Additional  damping  could  be  generated,  if  needed,  through  the  use  of  a 
large  dashpot  in  parallel  with  the  actuating  linkage.   However,  such  a 
dashpot  was  not  apparently  needed,  as  was  shown  through  investigation  of 
the  inherent  damping  of  the  rudder. 

C .   RUDDER 

The  description  of  the  moments  involved  in  the  rotation  of  a  rudder 
about  its  pivot,  while  moving  through  a  fluid,  would  in  general  include 
functions  involving  the  zero,  first,  and  second-order  derivatives  of  rudder 
angle  with  respect  to  time.   That  is,  if  a  moment  balance  were  written  for 
a  rotating  rudder,  the  equation  could  in  general  take  the  following  form: 


M  =  f^(e)  +  f2(e)  +  f^(e) 


(25J 


where  f  ,  f^,  and  f  would  be  functions  of  some  yet-to-be-determined  form. 

1    ^         o 

Further  consideration  of  the  natures  of  f,  and  f  revealed  that  expressions 

for  these  two  functions  were  present  in  the  baseline  simulation.   If  f^  = 


f  =0,  the  general  expression  is  reduced  to 


M  =  f^(e*)  •  (26) 


■^5 


and  f  is  recognized  as  being  of  the  form 

f^  =  10  (27) 

where  I  is  the  polar  moment  of  inertia  of  the  rudder  and  its  fluid  added 
mass.   If  f^=f2=0,  f^  is  recognized  as  the  sum  of  all  6  -  dependent 
parameters  in  the  rudder,  including  the  hydrodynamic  torque  caused  by  the 
rudder  lift  and  the  difference  between  the  location  of  the  rudder's  pivot 
and  its  center  of  pressure,  and  the  e  -  dependent  portion  of  the  system 
friction. 

The  function  f  ( e)  was  assumed  equal  to  zero  in  the  baseline  simulation; 
that  is ,  no  q  -   dependent  effects  were  included  in  the  rudder's  equation  of 
motion  (Eqn.  8).   The  validity  of  this  assumption  was  examined  in  the  final 
portion  of  the  investigation  of  load  and  linkage  damping. 

1.   Existence  of  a  Non-Zero  Function  f2(®) 

Initial  considerations  of  the  problem  involved  the  search  for  any 
indication  that  a  non-zero  function  fyi^)    even  existed.   Intuitively,  such 
a  relation  should  exist  in  some  magnitude,  at  least  under  a  condition  where 
the  rudder  was  subjected  to  a  constant  turning  rate  while  immersed  in  a 
stagnant,  viscous  fluid,  since  power  would  necessarily  be  expended  to  keep 
such  a  structure  moving  under  these  conditions,  thus  resulting  in  some  level 
of  torque  being  applied  to  the  rudder  shaft.   This  applied  torque,  then, 
would  by  Eqn.  (25)  be  equal  to  f2(e),  and  therefore  f2  could  not  be  identically 
zero. 

Another  indication  of  the  existence  of  a  non-zero  f2  was  provided 
by  Landweber  [17],  who  indicated  that  forces  on  lifting  surfaces  in  unsteady 
motion  were  due  to  a  combination  of  inertia  effects  and  circulation,  and 
that  a  transient  variation  did  exist  in  the  circulation  around  a  thin  air- 
foil impulsively  given  an  angle  of  attack  a,  causing  a  transient  behavior 
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of  the  lift  coefficient  of  such  a  surface.   Additionally,  Theodorsen  [18] 
studied  the  forces  on  a  foil  undergoing  harmonic  oscillations  in  connection 
with  studies  of  control  surface  flutter  and  aerodynamic  instability,  and 
found  that  the  actual  lift  coefficient  was  related  to  a  quasi-steady  lift 
coefficient  corresponding  to  the  instantaneous  angle  of  attack  of  the  foil, 
and  time.   Thus,  the  existence  of  a  non-zero  function  f^(0),  while  not 
confirmed  by  the  above  references,  was  at  least  shown  to  be  plausible.   The 
writer  therefore  attempted  to  locate  any  work  that  had  been  done  in  the 
analysis  of  the  forces  and  turning  moments  on  airfoils  or  rudders  subjected 
to  constant  turning  rate. 

2.   Damping  in  a  Rudder  Immersed  in  a  Stationary  Fluid 

No  references  were  found  describing  the  resisting  torque  developed 
by  a  rudder  subjected  to  a  constant  rate  of  rotation  while  immersed  in  a 
viscous  fluid.   Accordingly,  the  writer  attempted  to  approximate  the  drag 
which  could  be  expected,  using  elementary  concepts. 


Figure  15.   Rudder  Construction  Assumed  For 

Stationary  Rudder  Torque  Calculation 
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The  rudder  was  approximated  by  a  series  of  very  small  flat  plates, 
assumed  to  be  separated  from  each  other  so  that  the  flow  could  move  around 
the  plates  (Figure  15) .   For  a  flat  plate  placed  in  a  stream  flowing  perpen- 
dicular to  the  plane  of  the  plate,  the  drag  force  F„  can  be  expressed  as 


^D   =    ^d{'^)^  (28) 


where  C  in  this  case  is  the  drag  coefficient,  p  is  the  density,  A  is  the 
area  of  the  small  plate,  and  the  velocity  V  is  related  to  the  radius  r  and 
the  angular  rate  of  rotation  6  by  V  =  rO.   Thus,  for  a  differential  area 


dAp 


pr2e2 
^^D  =  ^d  I  ~T-  j  dAp 

If  the  area  is  assumed  to  be  equal  to  the  span  of  the  rudder  S  times  the 
incremental  distance  dr,dA^  =  Sdr,  then 

dFo  =  %    (  ^V-  )  Sdr 

The  incremental  torque  dTj^  caused  by  F^^  developed  at  a  distance  r  from  the 
pivot  point  of  the  rudder  would  be 

/  peM  .  3,  (29) 

d\  =  S  (  V  )  Sr  dr 

Thus,  for  the  portion  of  the  plate  from  r  =  0  to  r  =  r^, 

C  pSr  4 

T  =  -V^  e2  .  (30) 

R      8 

For  a  flat  plate  in  a  normal  stream  with  a  large  S/dr  ratio,  Cj^  is 
approximately  1.95  [19].   For  the  two  rudders  of  a  typical  submarine, 
the  following  physical  parameters  are  representative: 
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Parameter 


Value 


S 

r   rear  portion 


10  Feet 
7  Feet 
r   forward  portion    3  Feet 


P ,   sea  water 


64  Ibm/ff 


The  entire  9  -  dependent  torque  would  be  equal  to  the  sum  of  the  effects 
of  both  upper  and  lower  rudders,  including  forward  and  aft  portions: 


Tq  =  2  (Tr 


FORWARD 


^  ^Raft^  • 


(31) 


Substituting  in  the  previously  mentioned  values,  it  is  found  that 


To  = 


,  in. -lb. 


(32) 


when  6  is  expressed  in  degrees  per  second.   If  9  =  10  deg/sec,  T-        =   8800 
^  or  a>         >      Bmax 

in.-lbf.   If  an  asyinmetric  square  analog  component  and  a  potentiometer  were 
placed  in  the  simulator  as  in  Figure  16,  this  source  of  damping  could  be 
modeled. 


Figure  16.   Analog  Circuit  Diagram  for  Introduction  of  Tq 
The  appropriate  potentiometer  setting  expression  for  simulating  T^  would  be 


Ter 


^34 


1006  13 
m 


(33) 
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Substituting  numerical  values, 

a^^  «  .968  X  10"-^  .  (34) 

Thus,  the  initial  estimate  of  available  rudder  drag  at  zero  velocity 
indicated  negligible  amounts  of  damping  were  present  in  the  rudder 
turned  at  a  constant  rate  in  a  stationary  fluid.   The  writer  felt  that 
the  results  were  inconclusive,  however,  since  the  assumptions  which  were 
made  in  the  calculation  of  T  were  very  rough  indeed.   The  rudder  was 
not  formed  of  many  small  strips  with  spaces  between  them;  nor  was  the 
flow  as  simple  as  that  described  by  the  model.   Such  a  complicated  flow 
as  would  be  present  along  the  surface  of  a  large  flat  plate  rotating  in 
a  quiescent  fluid  required  a  depth  of  analysis  beyond  the  scope  of  this 
study. 

3.   Damping  in  a  Rudder  Immersed  in  a  Moving  Stream 

A  literature  search  revealed  that  an  analysis  had  been  made  of  the 
torque  behavior  of  a  turning  rudder  immersed  in  a  moving  stream.   Okada  [20j 
published  a  paper  dealing  with  the  effect,  which  was  subsequently  expanded 
by  Kennard  and  Leibowitz  [21]  in  their  studies  of  rudder  performance  as 
affected  by  ship  turning,  steady  change  of  rudder  angle,  and  propellor  race. 

Kennard  and  Leibowitz  [22]  noted  that  the  twisting  moment  observed  on 
a  rudder  shaft  while  the  rudder  is  being  turned  may  differ  considerably 
from  its  value  when  the  rudder  is  stationary.   They  went  on  to  explain  that 
this  effect  was  due  to  circulatory  motions  around  the  foil  caused  by 
viscosity  effects  of  the  water,  which  subsequently  generate  vortex  sheets 
in  the  water  downstream  from  the  rudder's  trailing  edge,  and  cause  moment- 
modifying  variations  in  the  pressure  distribution  over  the  rudder's  surface. 
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Variations  from  steady  state  twisting  moments  by  as  much  as  thirty  percent 
were  mentioned  as  being  observed  in  rudders  turning  at  constant  rates. 
The  mathematical  difficulties  associated  with  the  problem  were 
simplified  by  several  approximations.   The  rudder  was  replaced  by  the  flat 
plate,  shown  in  Figure  17,  of  uniform  chord  C  and  infinite  vertical  span. 
The  water  was  assumed  to  approach  the  plate  at  a  uniform  rate  and  at  an 
angle  of  attack  a;  its  velocity  was  denoted  by  U.   The  rudder  twisting 
moment  per  unit  span  was  identified  as  Mj,,  and  was  assumed  positive  when 
operating  in  a  direction  to  decrease  a.   The  pivot  point  of  the  plate  was 
identified  as  H,  and  was  located  at  a  distance  hC  from  the  leading  edge, 
h  being  assigned  some  numerical  value  between  0.25  and  0.50.   The  rudder 
turning  rate  was  assumed  positive  for  increasing  angles  of  a. 


Figure  17.   Flat  Plate  Approximation  to  a  Rudder 
Turning  at  Constant  Rate 

Mu,  the  twisting  moment  per  unit  span,  was  placed  in  nondimensional  form 
by  dividing  it  by  (1/2)  C^U^,  to  form  the  dimensionless  moment  coefficient 
C^    An  expression  for  Cj^  was  then  derived,  with  C^^^  being  represented  as 
the  sum  of  a  quasi-stationary  part  C^^q,  and  a  9  -  dependent  part,  AC^j^. 
Thus, 


^MH  "  ^MHO  ""  ^^MH 


(35) 
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where 

^MHO  "  ^^4  -h)  2  -na  (36) 


(|  -h)(l-2h)-(i-h)Al      (37) 


and 

C  da 
^   =U  dT 

and  A  is  a  function  of  il  and  a.   The  effect  of  A  is  to  increase  the  value 

of  ACj^j^  as  a  increases,  as  is  shown  in  Figure  18. 

Kennard  and  Leibowitz  reported  that  the  expression  for  C,   was  not 

MHO 

entirely  correct,  due  to  the  neglect  of  a  nonlinear  effect  in  its  develop- 
ment.  That  nonlinear  effect  concerns  the  aft  movement  of  the  center  of 
pressure  in  a  rudder  as  a  is  increased,  which  causes  a  change  in  the  sign 
of  Cj^Q  with  increased  a.   Also,  they  indicated  that  values  for  AC|^|,  obtained 
from  model  tests  suggested  that  the  calculated  values  of  ACj^  were  somewhat 
too  high,  although  the  scatter  of  the  experimental  data  in  the  tests  was 
considerable.   They  further  concluded  that  the  general  trends  resulting 
from  turning  real  rudders  in  water  were  modeled  in  their  theory,  although 
values  for  ^C..ii  were  somewhat  too  high,  such  as  a  calculated  increase  in 
one  case  where  h  =  0.3  of  roughly  48  percent,  where  in  practice  the  increase 
due  to  AC.«ii  ranged  from  15  to  40  percent.   It  was  suggested  that  tests  with 
larger  rudders  would  be  worthwhile,  and  that  refinement  of  the  analysis  was 
needed,  expecially  in  including  the  effects  of  finite  rudder  spans. 
4.   Application  to  the  Simulation  Problem 

The  analysis  just  described  seemed  to  be  applicable  to  the  simulation 
problem.   Although  the  expression  for  Cj^„^  was  apparently  unrealistic,  the 

development  ofAC   was  apparently  useful  in  predicting  the  effect  of  rudder 
MH 

turning  on  the  observed  moments  on  a  rudder,  although  calculated  values 
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seemed  to  be  somewhat  higher  than  values  which  were  experimentally  obtained. 
In  spite  of  its  inherent  shortcomings,  incorporation  of  this  modeling  of 
the  effect  of  rudder  turning  rate  on  the  observed  twisting  moment  seemed 
justified  in  light  of  the  reported  magnitudes  of  change  in  twisting  moment 
observed  in  practice. 

In  order  to  apply  the  above  analysis  to  the  simulation,  several 
assumptions  were  made.   It  was  assumed  that  a  was  identical  to  ^  ;  that  is, 
that  the  fluid  flowed  parallel  to  centerline  of  the  submarine.   Analysis  of 
drawings  of  the  SSN-671  submarine  rudder  revealed  that  h  was  approximately 
equal  to  0.28;  in  order  to  compensate  for  the  reported  high  nature  of  the 
calculated  AC„  ,  and  to  simplify  calculations,  h  was  assumed  to  be  1/4, 
thus  eliminating  dependence  on  A.   Also,  the  sign  convention  of  M„  was 
reversed  to  conform  to  the  simulation's  convention.   The  above  assumptions 
when  used  in  the  expression  for  AC   resulted  in 

Re-dimensionalizing  to  form  AM   the  change  in  twisting  moment  per  unit 
span  due  to  9,  and  then  multiplying  by  the  span  length  S,  r-esulted  in  the 
expression 

m^s  =  -  I  pc^sve  .  (39) 

For  a  SSN-671  rudder,  the  average  value  of  C  is  about  10  feet,  and  S  is 
also  about  10  feet.   In  a  typical  submarine  V  can  vary  from  zero  to  in 
excess  of  35  ft/sec,  and  0  was  assumed  in  the  baseline  simulation  to  vary 
from  zero  to  10  degrees/sec.   Assuming  p =  64  Ibm/ft^,  the  maximum  value  of 
the  twisting  moment,  (AM^^S),^^^  based  on  Vo  =  35  ft/sec,  would  be 


CAMoS)    =  5.72  X  10^  in-lbf 
^   rl  ■'max 


(40) 
max 


f4 


Typically,  a  submarine  has  two  rudders,  of  approximately  the  same  design, 
The  total  twisting  moment  available  would  therefore  be  twice  the  above 
value,  or  1,144  X  10   in-lbf.   An  expression  for  AM  S/(AM  S)    would  be 


AM^S 


(AMj^S) 


max 


-  (f  ]  (|-  )  Ml) 

\  o  /  »  ^max  / 


AMj^S/(  MjjS)^^^  is  dependent  only  on   for  a  given  value  of  V,  and  therefore 

could  be  simulated  as  damping  across  integrator  25  as  in  Figure  10.   The 

appropriate  expression  for  the  potentiometer  setting  a_„  would  be 

(  M,,S)    (V/V  ) 
H  max    o 
a   »  ; (42) 

le    (lOO)e 
max 

For  V  =  V  ,  and  the  conditions  otherwise  assumed  above,  since  V  appeared  to 

be  set  equal  to  V   for  the  other  V  -  dependent  parameters  in  the  baseline 

simulation,  a-_  =  0.1259.  A  level  of  0.12  —  feedback  was  previously  found 

m 
to  be  sufficient  to  approach  a  critically  damped  behavior  of  the  simulation, 

and  therefore  eliminate  the  oscillations  from  the  system.  Apparently,  then, 

a  theoretical  justification  for  the  inclusion  of  that  much  damping  in  the 

system  had  been  found.   The  fact  that  AC   was  dependent  on  V  did  not  allow 

Mn 

the  extension  of  the  above  damping  levels  to  lower  velocity,  since  as  V 
decreases,  so  must  the  level  of  damping.   The  authors  felt  that  the  damping 
effect  somehow  assumed  a  significant  value  even  at  low  or  zero  velocities, 
thus  indicating  the  existence  of  some  unexamined  or  undiscovered  contributing 
effects  such  as  those  previously  discussed,  which  would  be  independent  of  V. 
More  work  in  the  area  of  rudder  twisting  moment  and  lift  variance  with  rudder 
turning  angle  and  other  parameters  obviously  would  be  useful. 
5.   Summary 

The  investigation  of  linkage  and  load  damping  resulted  in  the 
location  of  a  significant  source  of  9  -  damping  in  the  motion  of  a  rudder 
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inunersed  in  a  moving  fluid.   A  significant  amount  of  damping  was  evidenced 
by  the  large  variation  of  the  observed  twisting  moment  produced  by  the 
turning  of  a  rudder  from  that  which  would  be  observed  if  the  rudder  were 
held  stationary  at  a  given  angle  of  attack.   On  the  basis  of  rudimentary 
analysis,  no  large  amounts  of  viscous  damping  were  located  in  the  packings, 
or  the  crosshead  structure,  although  a  slight  amount  of  damping  could  be 
caused  in  the  crosshead  through  redesign  of  the  structure.  A  very  rough 
estimate  of  the  twisting  moment  observed  on  a  rudder  immersed  in  a  stationary 
fluid  indicated  that  little  available  damping  existed  unless  the  rudder  is 
moved  through  the  fluid.   More  analytical  and  experimental  research  needs 
to  be  done  concerning  the  transient  lift  and  twisting  moment  characteristics 
of  moving  and  stationary  rudders  in  fluid  flow  fields. 
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IV.   RECOMMENDATIONS  FOR  FURTHER  WORK 

Several  topics  were  identified  as  requiring  investigation  during  the 
progress  of  this  project.   Due  to  time  limitations,  none  of  them  were 
pursued  in  detail.   Subjects  uncovered  in  the  initial  investigation  which 
remained  significant  at  the  end  of  the  period  that  this  report  covers 
included  the  need  for  a  comparison  of  the  actual  and  simulator  performances 
under  varying  conditions,  including  a  wide  range  of  ship  velocities. 
Studies  of  the  effect  of  using  a  higher  order  servovalve  and  uncompensated 
actuator  models  in  the  simulation  are  necessary.   Possibly,  the  re-checking 
of  the  consistency  of  the  numerical  values  used  in  the  baseline  simulation 
is  warranted,  although  the  discovery  of  a  significant  source  of  system 
damping  seemed  to  make  this  unnecessary.   Additionally,  the  investigation 
of  linkage  and  load  damping  emphasized  the  need  for  continuing  research 
into  unsteady  rudder  system  performance  characteristics,  including  rudder 
lift  and  twisting  moment  variations  with  such  variables  as  rudder  turning 
rate,  ship  velocity,  and  ship  turning  rate. 
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V.   CONCLUSIONS 

1.  Analysis  of  unsteady  rudder  behavior  identified  a  significant 
source  of  rate-dependent  torque  which  had  been  neglected  in  the  development 
of  a  submarine  steering  control  system's  analog  simulation.   Programming 

of  representative  levels  of  feedback  in  the  simulator  virtually  eliminated 
oscillations  in  rudder  rate  and  load  pressure  drop  which  were  thought  to 
be  unrealistic. 

2.  The  actual  performance  characteristics  of  a  submarine  steering 
control  system  under  varying  operational  conditions  appeared  to  be  loosely 
defined,  and  quantitative  comparisons  of  system  performance  seemed  to  be 
necessary  in  order  to  continue  the  development  of  a  realistic  simulator. 

3.  Studies  showed  that  the  inclusion  of  more  complex  models  of  the 
servovalve  and  actuator  could  have  important  effects  on  the  flexibility 
and  realism  of  the  simulator,  and  that  further  work  should  be  done  to 
accomplish  this  end. 

4.  Continuing  theoretical  and  experimental  research  in  the  field  of 
rudder  behavior  under  unsteady  conditions  was  suggested  in  order  to  improve 
the  current  concepts  employed  in  steering  system  design  and  simulation. 
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Appendix  A 
Explanation  of  Notation 


1.   Sat  (  X  ) 


uL 
LL 


Sat(x)^j 


ul   x>ul 
X   ll^x4ul 
11  x<ll 


This  can  be  represented  graphically  as  follows 


2.   Sgn(x) 


Sgn(x)  = 


+  1  x>0 

0  x=0 
-1   x<0 


This  can  be  represented  as  follows 


4-1 


-1 


69 


LIST  OF  REFERENCES 

1.  Morse,  Allen  C. ,  Electrohydraulic  Servomechanisms.  pp.  63-64.  McGraw- 
Hill,  1963.  

2.  Merritt,  H.  E.,  Hydraulic  Control  Systems,  p.  215,  Wiley,  1967. 

3.  Johnson,  G. ,  NAVSHIPRANDCEN,  Code  2762.3,  Personal  Communication. 

4.  Morse,  op.  cit.,  p.  142. 

5.  Johnson,  op.  cit. 

6.  Id. 

7.  Johnson,  Guy,  Unclassified  letter  to  R.  H.  Nunn  describing  the 
simulation,  20  December  1971. 

8.  Smith,  Leeland,  Code  260.4,  Mare  Island  Naval  Shipyard,  Vallejo,  Ca. 
94952,  Personal  communication. 

9.  Keller,  G.  R. ,  Hydraulic  System  Analysis,  pp.  139-142,  Industrial,  1969. 

10.  Kendall,  W.  T. ,  Code  260.41,  Mare  Island  Naval  Shipyard,  Vallejo,  Ca. 
94592,  Personal  Communication. 

11.  Chase,  Leonard,  Code  260.41,  Mare  Island  Naval  Shipyard,  Vallejo,  Ca. 
94592. 

12.  Merritt,  op.  cit.,  p.  34. 

13.  Ibid.,  p.  42. 

14.  Laughaar,  H.  L. ,  Steady  Flow  in  the  Transition  Length  of  a  Straight  Tube, 
Journal  of  Applied  Mechanics,  June  1942,  A55-A58. 

15.  Krieth,  F.,  Eisenstadt,  R.,  Pressure  Drop  and  Flow  Characteristics  of 
Short  Capillary  tubes  at  Low  Reynolds  Numbers,  ASME  Transactions,  July 
1957  1070-1078. 

16.  Shapiro,  H. ,  Siegel,  and  Kline,  J.,  Friction  Factor  in  the  Laminar  Entry 
Region  of  a  Smooth  Tube,  Proc.  2nd  U.S.  Natl.  Congress  of  Applied 
Mechanics,  June  1954   733-741. 

17.  Streeter,  V.,  Handbook,  of  Fluid  Dynamics,  1st  Ed.,  pp.  13-35,  13-36, 
McGraw-Hill,  1961. 

18.  Theodorson,  T.,  General  Theory  of  Aerodynamic  Instability  and  the 
Mechanism  of  Flutter,  NACA  Technical  Report  No.  496,  1949. 

19.  Hughes,  W. ,  Brighton,  J.,  Shaum's  Outline  of  Theory  and  Problems  of 
Fluid  Dynamics,  p.  85. 


70 


20.  Okada,  S.,  Investigation  of  the  Effect  of  the  Angular  Velocity  of 
Steering  Upon  the  Performance  of  Rudders,  Technical  Research 
Laboratory,  HITACHI  Shipbuilding  and  Engineering  Company,  Ltd.,  1958. 

21.  Kennard,  E, ,  Leibowitz,  R. ,  Theory  of  Forces  and  Moments  on  Rudders  as 
Affected  by  Ship  Turning,  Steady  Change  of  Rudder  Angle,  and  Propellor 
Race,  NSRDC  Acoustics  and  Vibrations  Laboratory  Technical  Note 

AVL-2  38-9A0,  pp.  47-71,  January  1969. 

22.  Ibid.,  p.  47. 


71 


INITIAL  DISTRIBUTION  LIST 


Officer-in-Charge 

Annapolis  Laboratory 

Naval  Ship  Research  and  Development  Center 

Annapolis,  Maryland  21402 

Attn:  Mr.  G.  L.  Johnson,  Code  2762.3  12 

Library 

Naval  Postgraduate  School 

Monterey,  California  93940  2 

Dean  of  Research  Administration 

Naval  Postgraduate  School 

Monterey,  California  93940  2 

Director 

Defense  Documentation  Center 

5010  Duke  Street 

Alexandria,  Virginia  22314  20 

Robert  H.  Nunn,  Code  59Nn 
Naval  Postgraduate  School 
Monterey,  California  93940  5 

Mr.  W.  J.  Blumberg,  Head 

Simulation  and  Analysis  Branch 

Code  A622 

Naval  Ship  Research  and  Development  Center 

Annapolis,  Maryland  21402  1 

Commanding  Officer 

Navy  Nuclear  Power  School 

Mare  Island,  Vallejo,  Ca.   94592 

Attn:   Ensign  Ernest  E.  Wessman  1 

Mr.  Ralph  Leibowitz,  Code  194 

Naval  Ship  Research  and  Development  Center 

Bethesda,  Maryland  20034  1 

Mr.  William  T.  Kendall,  Code  260.41 

Mare  Island  Naval  Shipyard 

Vallejo,  Ca.   94592  1 

Professor  G.  Thaler,  Code  52Tr 

Naval  Postgraduate  School 

Monterey,  California  93940  1 


72 


Security  Clonlficttlon 


DOCUMENT  COHTROL  DATA  -R&D 


(Stcurliy  cltnlllcailon  of  (/(/•,   body  ol  mbtlraci  tind  IndtMing  tnnolallon  mu*l  bt  »nl»r»d  wh»n  lh»  ovrtll  report  li  cl*t»t(le4) 
I     OWiOiNATiNC   ACTIVITY  (Corporaf  author)  la,  REPORT  i«CU»<lT¥    C  L  Al«i  r  ic  A  TION 


1 


Naval  Postgraduate  School 
Monterey,  California  93940 


Unclassified 


lb.    OROUP 


»      RtPOPI  T     Tl  TL  1 


Analysis  and  Simulation  of  a  Submarine  Steering  Control  System 


4.   DKtCRiPTivK  NOTES  (Typ»  ol  rapori  and,lnelu»lv»  daft) 


I     Au  THORiSi  (Flral  rtama,  mlci<Jla  Initial,  laat  name; 

Ernest  E.  Wessman 
Robert  H,  Nunn 


(.  WCPORT  DATE 

1  August  1972 


Ta.    TOTAL  NO.  OF  PAOC* 

74 « 


7b.  NO.  or  NKFt 

22 


•  a.    CONTRACT    OR    ORANT    NO. 


6.    ^MOJICT  NO. 


ft«.    ORIQINATOR'S    REPORT   NUMBER!*) 


NPS-59NN72081A 


Annapolis  Laboratory 
Naval  Ship  Research  and  Development 
Center,  WR  2-0712 


»6.  OTHER  REPORT  NO(*l  (Any  othor  numbata  that  ansy  ba  aaalpyad 
thia  raport) 


10.    OltTRieuTION    STATEMENT 


■Diotrib-ution  Ij 

Proprietary  information,  8  August  1972. — Other  requests- 

for  thio  dooumcnt  must  b-e  ^i°-&g^rred  to  Naval  Postgraduate 
Q^^ee4-^_M^-il:crcv.  Callfuiuid  939A(>T-&od-e-e^-. 


n.    SUPPLEMENT ARV    NOTES 


DOWNGRADED 
APPROVED  FOR  PUBLIC  RELEASE 


12.    SPONSORING   MILl  TAR  Y    ACTIVITY 


Annapolis  Laboratory 
Naval  Ship  Research  and  Development 
Center 


13.  AKSTRACT 


An  existing  computer  simulation  of  a  hydraulic  submarine  steering  control  system 
was  analyzed  to  determine  the  source  of  its  unrealistic  behavior,  which  included 
verifying  mathematical  equations  and  numerical  parameters,  duplicating  the  simula- 
tion on  another  computer,  and  investigating  the  validity  of  the  assumptions  made 
in  its  development.   Rudder  turning  rate  and  load  pressure  drop  oscillations  which 
were  thought  to  be  unrealistic  were  caused  by  the  neglect  of  significant  rudder 
moments  related  to  the  control  surface  turning  rate.   Knowledge  of  actual  system 
performance  was  found  to  be  incomplete  and  further  research  was  suggested.   Poten- 
tially beneficial  modifications  in  the  servovalve  and  actuator  simulators  were 
identified.   Continuing  work  in  the  field  of  unsteady  rudder  behavior  was  suggested 
in  order  to  improve  simulator  and  actual  system  designs. 


l#ll^    <  M9V  tt 

t/K  oiei-607-oin 


1473       (PAG£   1) 


73 


"i^iliy  O&talHM'tl' 


M 


A-»l«9» 


SeLuritv  Classification 


KEY     AORD5 


RO  L  E  *n 


L  {  N  K      C 


Steering 

Simulation 

Hydraulic 

Rudder 

Submarine 

Control 


DD  ,Z':,AA73   BACK, 


'N      .■  1  0  1  ■  80  7  -  68  ?  1 


74 


Security  Classification 


Al 


DUDLEY  KNOX  LIBRARY  -  RESEARCH  REPORTS 


5  6853  01068097  8 


Shelved  numerically  in  tlie  vav 


